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Abstract 

We consider the problem of global stability of solutions to a class of semilinear wave equations 
with null condition in Minkowski space. We give sufficient conditions on the given solution which 
guarantees stability. Our stability result can be reduced to a small data global existence result for 
a class of semilinear wave equations with linear terms B M "9 M $(t, x)d v (j>, L M (t,x)c' M and quadratic 
terms h^"(t,x)d l j,(j>dv(f) where the functions <fr(t,x), L M (i,a;), W^t^x) decay rather weakly and the 
constants B^ v satisfy the null condition. We show the small data global existence result by using 
the new approach developed by M. Dafermos and I. Rodnianski. In particular, we prove the global 
stability result under weaker assumptions than those imposed by S. Alinhac. 

1 Introduction 

In this paper, we study the behavior of solutions to the Cauchy problem 
ju w =J\f(dw), 

\w(0,x) = $ (x) + e</> (x), d t w(0,x) = $i(x) + e<j> 1 (x) 



(1) 



in Minkowski space with initial data $j(x), 4>i(x) € Cq°(R 3 ). The nonlinearity Af(dw) is assumed to 
satisfy the null condition, that is, Af(0) = DAf(0) — and the quadratic part of Af(dw) is A a @ d a wdpw 
with constant coefficients A a @ such that A a P£ a £p = whenever £jj = £f + £f + Cl- 
in [1], S. Alinhac studied the stability of large solutions to the quasilinear wave equations 

j a w + g a ^d 7 w ■ d afs w = 0, 
\w(0,x) = $o(x), d t w(Q,x) = $ 1 (x) 

in Minkowski space, where g al3 ^ are constants satisfying the null condition ( see [11]). More specifically, 
starting with a global solution x) G C°°(R 3+1 ), consider the Cauchy problem with perturbed initial 
data ($(0, x) + e<fio, dt$(0, x) + e</>i). He showed that if $ satisfies the condition 

<r r '>-'l'<^, • c 2 . 5^ |r*a*| < <7 (1 + + ]r - (2) 

i=l |fe|<7 

for some positive constants ccq < 1 and Co, then the solution exists globally and is close to $. Here 
r denotes the collection of Lorentz vector fields, see [10]. The problem of global stability of $ can be 
reduced to the following small data Cauchy problem 

f n<f> + g a ^d 7 cj> ■ d aP cf + g^'d^dp^ + g^-rdfrGdat = 0, 
|</>(0,x) = e0 o (x), d t <f>{0,x) = e0i(x) 

with given function $ satisfying condition (2). The approach in [1] relies on the vector field method. In 
particular, S. Alinhac used the scaling vector field S = tdt+rd r with weights growing in t as commutators. 
The use of such weighted vector fields requires one to make the rather strong assumption that the given 
solution $(t,x) decays uniformly in time t as in (2). 
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In this paper, we use the approach developed in [5], [22] to treat the problem of global stability 
of solutions to nonlinear wave equations, we use a new method for proving decay for linear problem, 
developed by M. Dafermos and I. Rodnianski in [5]. This new method avoids the use of vector fields 
containing positive weights in t, e.g., S — td t + rd ri Li = Xid t + tdi. Traditionally, the vector fields from 
the set r, including S, Li, are used as multipliers or commutators. The new approach only commutes the 
equation with d t , Xidj — xjdi and allows us to obtain the stability results under conditions on $ weaker 
than those imposed by inequalities (2). We now describe the assumptions and the main results. 

We denote (d t ,d Xl ,d X2 ,d X3 ) by d and (d Xl , d X2 , d X3 ) by V under the coordinates (t,x\,X2,Xs). We 
also use the null coordinates u = v = defined by the standard polar coordinates (t, r, co) in 
Minkowski space. The vector fields, used as commutators , are 

Z {H^j, — ^ij — Xidj Xjdi, 

where Greek indices run from to 3 while the Latin indices run from 1 to 3. 
Definition 1. We call $ e C°°(R 3+1 ) a {6,a,t , R 1 ,C )-weak wave if 

(i): |0$(t,s)| <C Q , t<t Q , 

(ii): \d$(t,x)\ <C (l + r)-i(l + (t-\x\)+)-?- 4a , \x\>R u t>t , 

(Hi): \d v $(t,x)\ < C (l + r)- 1 - 3a , \x\ > R u t>t , 

(iv): \d${t,x)\ < 5a{\ + r)- 1 -°, \x\ < R u t>t 

for some positive constants 5, a, to, Ri, Co, where d v — d t + d r . Here (t — \x\) + — max{0, t — \x\}. Without 
loss of generality, we assume a < j and R\ <t . 

Remark 1. Solution of a free wave equation in Minkowski space □$ = with compactly supported 
initial data decays uniformly in time t and is always a (5, a, to, Ri, Co)-weak wave for some constants 
S,a,to,R\,Co . We remark here that a weak wave does not have to decay uniformly in time t in the 
cylinder {(t, x)\\x\ < Ri}. 

In our argument, we estimate the decay of the solution with respect to the foliation E T , defined as 

S T := {u = u T ,v > v T }, 
S T := {t = r,r < R}US T , 

where u T = v T = The radius R is a to-be- fixed constant. The corresponding energy flux is 

E[4>](t):= [ \d(t>\ 2 dx+ ( (\d v (f>\ 2 + \f^\ 2 ) r 2 dvduj, 

Jr<R JS T 

where J/ denotes the induced covariant derivative on the sphere of radius r. We denote 

Eo= J2 [ \dZ k 4>(0,x)\ 2 dx, 

where 0(0, x) = <fio{x), dt<f)(0,x) — (f>i(x). Here k stands for multiple indices, namely if k = (feijfe), then 
Z k = Qfei T fe2 5 q = And if k < k, then k x < k u k 2 < h- 

In addition to the assumption that the nonlinearity J\f(dw) satisfies the null condition, we assume M 
is smooth and 

J\f(d<& + d<t>) = W(<9<I>) + A^d^dut + Af»(d$)d li 4, + N^{d<5>)d^d v 4> + 0(|^| 3 ), 
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when d(p is small. The coefficients J\f»(d$), Af^(d$) satisfy 

\ZPAf»(d$)\ < C(d$) \Z' 3 ' d<5>\ 2+ao , V|/3|<4, 

l/3'l<l/3| (3) 

\z Ar» v (d$)\ < c(9$) ^ |z"'0*| ao , v|^|<4 

l/3'l<l/3| 

for some positive constant ao- The constant C(<9$) depends only on J2 1 1 -^c^ lie - 

\p\<4 

We now state our main results. 

Theorem 1. Suppose the nonlinearity N{dw) satisfies the null condition and condition (3). Let $ € 
C°°(IR 3+1 ) be a solution of (1) whene~0. Assume Z k <& is (5, a, to, R\,Co)-weak wave, V|fc| < 4. Suppose 
the initial data 4>o(x), 4>\{x) are smooth and supported in {\x\ < Rq}. Then there exists 6q > 0, depending 
only on the constants A a ^ , and e > 0, depending on E , R , A a ^ , cto, a, t , R\, Co, such that for any 
5 < So, e < eo, there exists a unique global smooth solution w of equation (1) with the property that 3 
positive constant R, depending on to, ex., ao, Ri, Co, Ro, such that for the foliation S T with radius R, 
the difference (f> — w — $ satisfies 

(1) Energy decay 

E[4>](t) < CEoe 2 (l + T)- 1 -? a ' , a' = min{^,a}. 

6 

(2) Pointwise decay: 

\4>\ < C^/Eo~e(l + r)-\ 
]T \d^\<C^E~oe(l + r)-^(l + \t-r + R\)-^-^, a '-min{^,a}, 

l/3|<2 

where C depends on R, ao, a, to, R\, Co- 

Remark 2. The weak decay of d$> in the spatial direction ( (1 + |^|) ^ ) excludes general cubic nonlin- 
earities of N{d<&) ( cubic nonlinearities satisfying the null condition are allowed). However if condition 
(ii) in the definition of weak wave $ is improved to 

\dZ $\ < C (l + r)-3" a (l + (t - \x\) + )~^- ia , 

l/3|<4 

then it is sufficient to assume 

\Z f> Afi l (d$)\<C(d$) \ z<i ' d ^?i V|/3|<4. 

I0'l<|j8| 

This allows any cubic( or higher) nonlinearity of Af(dw). 

Since $(f , x) solves (1) for e = 0, the problem of global stability of <I> is then reduced to the following 
small data Cauchy problem 

{u<l> + N(*,<i>) + L{d<l>) = F{d<l>), 
\<f)(0,x) = e<j> (x), <j> t (0,x) = e<j> 1 (x), 

where N (<&,</>) = B a ^ Q $ • dp<f>, L(d<f>) = L' i (t,x)d tl <j). The nonlinearity F{d<j>) is of the form 

F(d<p) — A^ v d^<pd v (j) + Q(d(j>) + cubic and higher order terms of def), 
Q{d(j>) = h^(t,x)d^d v (t>. 
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Here A^ v , B^ v are constants satisfying the null condition [11] and $(t, x), L^(t,x), h^ v (t,x) are given 
functions. For the stability problem, we have = —2A^ V and <j> = w — <&. However, it is of inde- 
pendent interest to consider the above small data Cauchy problem with linear terms N($>,<fr), L^d^cj) 
and quadratic terms h^ v (t, x)d^cf)d l/ 4 l where the functions $(t, x), L^(t, x), h^ v (t, x) decay rather weakly, 
given as follows: 

For positive constants 6, a, to, Ri, Co, we assume Z k <& is (5, a, t 0} R\, Co)-weak wave, V|fe| < 4 and 

\d 2 ZP$\ < Co, V|/3| < 2. 

Similarly, we assume 

\dZ p L^\ < C , V|/3| < 2. 
For t < t , we assume Z^L^it.x), Z^h^ v (t,x) are bounded, that is, 

iZ^ftaOI + IZ^frx)! <C 0) V|/3|<4, Vt < t . 

For t > t j w e assume 

|Z^(t,x)| <C (l + M)-i", V|/3| < 4, Vt>t 
and L^(t,x) satisfies one of the following two conditions 

\Z?L> 1 {t,x)\<5a(l + \x\)- 1 -Z a , V|/3|<4, Vt > i (5) 

or 

|Z"L"(t,aO| <C (l + |.T|)- 1 - 3tt (l + (t-| a; |) + )- a , V|/3|<4, Vt>t . (6) 
Theorem 1 follows from: 

Theorem 2. Let <$>(t,x), L^(t,x), h^ v (t,x) be given smooth functions satisfying the above conditions. 
A^ v , B^ v are constants satisfying the null condition. Assume the initial data 4>q{x), <j)\(x) are smooth 
and supported in {\x\ < Rq}. Then there exists <5 > 0, depending only on the constants B^ v , and e > 0, 
depending on Eq,R$, A^ v ' , B^ v , a, to, R\,Cq, such that for any S < Sq, e < eo, there exists a unique 
global smooth solution <f> of the equation (4) with the property that 3 positive constant R, depending on 
to,a,Ri,Co,Ro, such that for the foliation S r with radius R, the solution <f) satisfies 

(1) Energy decay 

E[$\{T)<CEoe 2 {l + T)- 1 - 1 * a . 

(2) Pointwise decay 

H<Cv/^e(l + r)-\ 
J2 ^ CV^e(l + r)-*(l + |t-r + iJ|)-*-*, 

|/3|<2 

where C depends on R, ao, a, to,Ri, Co- 

Remark 3. Notice that a can be arbitrarily small. The decay assumptions on L^(t,x)( condition (5)J 
and W(t,x) are sharp in the sense that there exists soliton solution( independent of time t) to the linear 
wave equation if L^(t,x) behaves like (1 + |a;|) _1 and any nontrivial C 3 solution of the equation 

with compactly supported initial data blows up in finite time [9]. 
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Remark 4. We can also consider equation (4) with zeroth order linear term L n (t,x)4>, leading to the 
same conclusion provided that Lo(t,x) decays like (1 + |a;|)~ 3 ~ a . Hence for the stability problem of large 
solution( Theorem 1), specific dependence on w of the nonlinearity M(w) is also allowed. 

Remark 5. For simplicity, we consider the equations in Minkowski space. However, as in [22], the same 
conclusion holds on curved background (R 3+1 ,g) with metric g merely C 1 close to the Minkowski metric 
and coinciding with the Minkowski metric outside the cylinder {(t, x)\\x\ < R}. 

Remark 6. It is not necessary to require that the initial data have compact support. The general as- 
sumption on the initial data can be that the following quantity 



The small data global existence result of nonlinear wave equations satisfying the null condition in 
Minkowski space was first obtained by D. Christodoulou [2] and S. Klainerman [11]. The approach of [11] 
used the vector field method, introduced by S. Klainerman in [10]. Various applications of the vector field 
method to nonlinear wave equations could be found in [12], [13], [15], [17], [19], [20] . In particular, the 
celebrated global nonlinear stability of Minkowski space has been proven by Christodoulou-Klainerman 
[3] and later by Lindblad-Rodnianski [14]. 

The main difficulty of considering nonlinear wave equation (4) with linear terms B tJ,ly d^(t,x)d u (f>, 
L^(t,x)d ll (j) an d quadratic terms h^ v (t, x)d ll (j)d„(f) is the rather weak decay of the functions $(t,x), 
L^(t,x), W(t,x). Previous works have relied on the fact that these functions decay to zero uniformly 
in time t, which is not necessary in this context. In fact, we even allow these functions to stay static ( 
independent oft) in the cylinder {(i, x)||x| < R\}. Although we require S to be sufficiently small, which 
is the only smallness assumption here, 5 depends only on the constants B^ v . 

Our argument here is similar to that in [22], which relies on a new approach, developed by M. 
Dafermos and I. Rodnianski in [5] . This new approach is a combination of an integrated local energy 
inequality and a p-weighted energy inequality in a neighborhood of the null infinity. However, due to the 
weak decay of <9$, we are not able to obtain the integrated local energy inequality and the p-weighted 
energy inequalities separately as in [22]. We thus consider these two inequalities together, see Proposition 

1 in Section 2 for details. 

The plan of this paper is as follows: we establish an integrated energy inequality in the whole space 
time and two p-weighted energy inequalities in Section 2. In Section 3, we use Proposition 1 to obtain the 
decay of the energy as well as the pointwise decay of the solution under appropriate boostrap assumptions; 
in the last two sections, we close our boostrap argument and conclude our main theorems. 

Acknowlegements The author is deeply indebted his advisor Igor Rodnianski for suggesting this 
problem. He thanks Igor Rodnianski for sharing numerous valuable thoughts as well as many helpful 
comments on the manuscript. 

2 Notations and Preliminaries 

In Minkowski space, we recall the energy-momentum tensor 




is finite. 



= dfi<j)d u (f) - -g^d^^djej). 



Given a vector field X, we define the currents 



K x [<t>]=T^[<P}n l 



.x 
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where -k x v = is the deformation tensor of the vector field X. Recall that 

D» J* [<j>]=X(ci>)n<f> + K x [<{>}. 

Let n be the unit normal vector field to hypcrsurfaces, da the induced surface measure. We denote (ivol 
as the volume form in Minkowski space. In null coordinates, we define the null infinity from t± to T2 as 
follows 

T T \ := {(u,v,u)\u n <u< u T2 ,v = 00}. 
The corresponding energy flux is 



which is interpreted as a limit when v — > 00. Denote 

E[<j>]{T) = E[<p]{T)+mi. 

Here E[<P](t) is defined in the introduction with the foliation E T . 
Taking a vector field 

X = f(r)8 r , 

for some function f(r), consider the region bounded by the hypersurfaces S Tl and E T2 . By Stoke's 
formula, we have the identity 

f J*[<p}n»da- f J*[cj>\n»dcj- [ ^J*[cj>\n»dcj 
= [ 2 l D»J*W™\= f 2 [ nd>.X(cl>)+K x l<f>}dvol, (7) 



where 



K X [<t>] = TH4>]<„ = (\f + r-'DO^r + (\f r-\f)(d r( f>f \fW^ 



We use /' to denote d r f throughout this paper. 

Choose another function \ of r, we have the equality 

- X d^d l( j> + in x • <f> 2 = Id" (d^ x ■ 2 - xd^ 2 ) + x4>nct>. 

Add the above equality to both sides of (7). Define the current 



Then we get 



f jX[4>]n»da- I J x [<f>}n^da- lj x \^do 

= fj £ (fd r <P + <Pxp(t> + (r- 1 / + \f - x)(m 2 

+ (X r-\f + \f)(dA) 2 + ( X - \f)\n? ' 



(8) 
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(9) 



For X = T = dt in (7), we have the energy inequality 

E[<f>]fo) < E[<f>]{n) + 2 r f \n<p\\d t (j>\dvol. 

We now fix the radius R of the foliation £ T . First under the assumptions in Theorem 2, we choose new 
positive constants a', t' , R[ such that 

\d v Z k <5>\<5a'{l + \x\)- 1 - :ia \ t>t' , \x\>R[, \k\ < 4, (10) 

\0Z k <P\ < 25a' (1 + \x\y 1 - a ' , t>t' , \x\<R[, \k\ < 4, (11) 

(l + t' ) a 'Sa' >C . (12) 

In face, since Z fe $ is (S, a, to, R\, Co)-weak wave V|fc| < 4, choose a' — |a and R[ large enough such 
that (10) holds. Then (11) and (12) are satisfied if t' is sufficiently large. The other conditions are also 
satisfied for these new constants a' , t' , R[. Then the radius R can be fixed as 

R = t' Q + i?o , 

where Rq is the radius of the support of the initial data. To avoid too many constants, we still use the 
constants a, R±, to to denote a', R[, t' respectively in the sequel. 

The following several lemmas, which have been proven in [22], will be used later on. 

Lemma 1. On S T UX r 



' / \<f>\* 



duj < E[4>](t). 



Lemma 2. If (p is smooth, then 



Corollary 1. In the exterior region r > R 

f \d v (rcj>)\ 2 dvdu - f <& 
Lemma 3. Suppose f and \ satisfy 



\f\<C u \X\< 



1 + r' 



r dvduj 



\x'\< 



<2E[</>](t). 



C x 



(1 + r) 2 



for some constant C\, then 



< QCiE[(p](r). 



Remark 7. If E[4>](t) is finite, all the above statements are also valid if we replace E[tp](r) with (r). 
Finally, we denote 



(dct>, 



1 + r 



g(p,r) := / r p \d v ip\ 2 dvdij, g(p,r) := / r p \d v ip\ 2 dvdu>, 



T2 .. 
Tl ■ 



(1 + r) ^<?(p,r)dr, 



(1+r) Pg(p,T)dT, 



DP[F] T T l:= [ 2 f |F| 2 (l + r)^ +1 dvol 

Jt, is T 



for Vp > 0, (i € R 1 . Here ip = r<f>, d v — d t + d r and 8<p is a four dimensional vector with norm 
\d(j)\ 2 — \d<p\ 2 + ri+rf i similarly for 9^0. Throughout this paper, we use the notation A < B for the 
inequality A < C-B with some constant C, depending on R, A^ v ', B^ v , a, to, Ri, Co- 
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3 Weighted Energy Inequalities 

In this section, we use the multiplier method to establish an integrated local energy inequality and two p- 
weighted energy inequalities. The integrated local energy inequality was first proven by C. S. Morawetz 
in [18]. We follow the method in [4] to obtain the integrated local energy inequality here. In [5], M. 
Dafermos and I. Rodnianski introduced the p-weighted energy inequalities in a neighborhood of null 
infinity. These two estimates, which the decay of the energy flux E[4>]{t) relies on, were shown separately 
in [22]. Due to the weak decay of the functions $>(t,x), L^(t,x), we are not able to show these two 
estimates separately. We hence consider them together. 

Consider the following linear wave equation 

{n<f> + N(*,<t>) + L(d<i>) = F, 

[4>(0,x) = e<t>i(x), d t 4>(0,x) = e0i (x), 

where iV($, (/>) = B^d^ ■ d v (f), L(d<j>) = L»(t, x)d li (f>. We have the following key estimates. 

Proposition 1. Suppose that $ is a (5, a, to , Ri, Co) -weak wave for positive constants 5, a, to, Ri, Co- 
Assume the given functions L^(t,x) satisfy 

\L»(t,x)\<C , t<t 

and one of the following two conditions 

\L»(t,x)\ < fo(l + Izl)- 1 " 3 ", i>t 0) 

or 

\If(t,x)\ < C (l + |^|)-i-3« (1 + {t _ I^J )_,_)— «, vt > t Q . 

Suppose the constants a, to, 5, Co obey the relation (12). Then there exists <5 > 0, depending only on the 
constants B^ v , such that for all 5 < So, solution (p of equation (13) has the following properties: 

(1) Integrated local energy estimate 

T 2 r l^|2 



+ (1 + n)- 1 - 2 " (g(l + 2a, r 2 ) + jP" Ti a D 2a [F}?dT^j 



(14) 



(15) 



(2) Energy bound 

EMr 2 ) < EMn) +D a [F]% + (l + r^D^FY^ 

+ (l+ri)- 1 - 2a (g(l + 2a,T 2 ) + J\? a D 2a [F}?dT^j . 

(3) p-weighted energy inequalities in a neighborhood of null infinity 

9(1,t 2 ) + f * E[4>]{r)dr < g(l,Tx) + T^ a D 2a [F] T T \ + r{~ a E[<j>]{jx) + f '(1 + r)- a D 2a [F] 
J n J ri 

+ (l + r!)- 2Q r T L i a D 2a [F]?dT+(l + T 1 )- 2a g (l + 2a,T 1 ) 1 (16) 

Jti 

g(l + 2a, t 2 ) + G[0, 2a]£ < g(l + 2a, n ) + (1 + n) 1 -" E^n) 

+ H T? a D 2a [F]?dT + (l + n) 1+ i a D 2a [F} r T l. (17) 
J n 
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Remark 8. We mention here that variants and generalizations of estimate (14) can also be found in 
[16], [17]. 

The following corollary will be used to derive the energy decay estimates when commuting the equation 
with the vector fields Z. 

Corollary 2. Assume the given functions <fr(t,x), L^{t,x) and the constant S satisfy the conditions in 
the above proposition. Then for solution (f> of (13), we have estimates for N = B^d^ ■ d v (\), L — L^d^cj) 

D 2a [N]H + D 2a [L] T T 2 i < E[4>]{ T1 ) + D a [F]^ + (1 + Tl )- a D 2a [F]^ 

+ (1 + t 1 )- 1 - 2q (g(l + 2a,T 2 ) + T^ a D 2a [F] T /dT^j . 

Under appropriate boostrap assumptions on the nonlinearity F, the above inequalities lead to decay 
of the energy flux (t). We discuss the integrated local energy inequality and the p- weighted energy 
inequalities separately. And then combine them together to prove the above proposition. The following 
two lemmas will be used frequently. First define 

A = 10sup{|BH}. 

Lemma 4. Let N = B^d^d v <t>. Then 

\rN\ < A + \d v $>\\di[\ + \d$\\4>\), ^ = r4>. 

Proof. By our notations 

rN = rB^d^ ■ d v <j> = B^d^ ■ d„i[ - B^d^d v r ■ <j). 
The lemma then follows from the fact that B^ satisfies the null condition and the inequality \dr\ < 1. □ 

Lemma 5 (Gronwall's Inequality). Suppose A(t), E(t) are nonnegative functions on [t\,T2[ . Assume 
that E(t) is nondecreasing on this interval and (3 is a positive number. If 

A(t)<E[t]+C [ (l + s)- 1 - f! A(s)d s , VTe[T l7 T 2 ], 

J Tl 

then 

A^^exp^-^l + n)-^)^), Vre [r l7 r 2 ]. 
Proof. See [21]. □ 

3.1 Integrated Local Energy Inequality 

We follow the idea used in [4] by choosing appropriate functions / and x such that the coefficients on 
the right hand side of (8) are positive. The left hand side can be controlled by the energy flux E[<j)] by 
Lemma 3. We thus end up with an integrated energy inequality in the whole space time. We now discuss 
this in detail. 
Take 

/ = P - r, , w , X = r~ 1 f, P=-. 

(1 + r) a a 

Notice that 

(1 + r) a - 1 ^ a 
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We have 



r- 1 f+\f-X = X-r- 1 f+\f = j i ' 



(l + r) 1 +« : 



1 /3((l + r)"-l) _ 1 > 1 

^ O*' »h i / 1 i ~.\ 1 -4-rv — 



2 r(l + r) a (l + r) 1+Q ~ (1 + r) Q+1 ' 

1 a + 1 



2 A r(l + r) 2 +«' 
Hen the energy inequalities (8), (9) together with Lemma 3 imply that 

£ X (1 +%+« dY ° l ~ 12 ^Mti) + 13/?^" £ |F-Ar-Lp^vol. (18) 

To proceed, we have to estimate the linear terms AT($, 0), L(d<f)). We first consider the case t 2 > t\ > t . 
For t >t , notice that on E T 

C (l + + (* - M)+)~ Q < C (l + r)" a < C (l + to) - " < 5a 

by the inequality (12)(we have assumed this inequality in Proposition 1). Hence under the conditions on 
the functions L^(t,x) in Proposition 1, we always have 

f \L(d<p)\\dc/)\dvol < Sa [ * [ (1 + r) _1_a |00| 2 dvol. 

For 7V($,</>), we consider it inside and outside the cylinder {\x\ < Ri} separately. When r = \x\ < R\, 
the null structure of iV($, 4>) is not necessary. We has to rely on the smallness of S. Since $ is a weak 
wave, condition (iv) of Definition 1 implies that 

\N\ = \B al3 d a $d (t>\ < A5a{\ + r)- l - a \d(f)\. 

For r > Ri, the null structure of AT($, (f>) is of particular importance. By Lemma 4, it suffices to estimate 
the three terms r _1 |d<I>||<9„V||d0|, r _1 |9 t ,$||9r/'ll^|) r_1 |^H0ll^l- Without loss of generality, assume 
Ri > 1. For the second term, inequality (10) shows that 

Id^Hr-^Vll^l < 2Sa(l + r)- 1 - a \d(/)\ 2 . 

On S T n {\x\ > R\ > 1}, for the first term, we have 

r-^IMI^I < C{l + r)-'i{l + T)-^- Aa \d' v tP\\d(t)\ 

< 5a(l + r)- 1 -"^ 2 + C(l + T)- 1 - Sa r- 2+a \l%1>\ 2 . 

Here we denote C as a constant depending on a, R = to + Ro, A a/3 , B al3 , Co, 5. Similarly for the third 
term, we have 

r-^d&Mldtl < 8a(l + r)- l - a \B(t>\ 2 + C(l + T)- 1 - 8a r- 2+a \<l>\ 2 . 

It remains to control r~ 2+a \(f>\ 2 . We use the Hardy's inequality outside the cylinder {\x\ < R}. By 
Lemma 1, we have 



[ \i)\ 2 (T,v,w)dw < C [ \xp\ 2 (T,v T ,uj)duj + C ( I { \d v ip\doodv 

<CE[cj>]{T)+C f ( r 1+2a \d v ijj\ 2 dujdv f V r-^dv (19) 
< CE[4>](t) + Cg(l + 2a, r), (r, v, u) e S T , 
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where v T — -^J 1 , v = L ^-. Hence for all p < 1 + 2a 

f r p - 3a (f> 2 dvduj= f r p - 2 - 3a f \i(;\ 2 dujdv <CE[4>}(t) +Cg(l + 2a,T). 

J St <J V t J LO 

On the other hand, Lemma 2 shows that 

/ tfdvduj < CE[<f>]{r). 
Interpolate with (20) for p = 1 + 2a. We derive 

r a <p 2 dvdu < CE[<P](t) 1 -'* (e[4>](t) + . 9 (1 + 2a, t)) 7 
<CE[4>]{T)+CEmT) 1 -^g{l + 2a 1 T)\ 



(20) 



where 7 = jz^- This gives estimates for (f> 2 outside the cylinder < R}. 
In the region R\ <r <R, using Sobolev embedding and Lemma 1, we get 

f 4> 2 du<C { <f> 2 du +C f \d r (j)\ 2 dxleqCE[(j)}(T). 

Jui J cj r—R Jr<R 

Therefore we can estimate r~ 2+a \cf)\ 2 outside the cylinder {|x| < R\} as follows 

/ r- 2+a \(t)\ 2 da = [ r- 2+a (t> 2 dx+ f r a ^> 2 dvduj 

J {r>fii}nS T JR!<r<R J S T 

< CE[4>](t) + C£M(t)^j(1 + 2a, t)~< . 

Inside the cylinder {|x| < we use the assumption that <9<I> is small. Hence combining the above 

estimates, we can bound the linear term (\N(^,<p)\ + \L(d(j))\)\d(j)\ in (18) as follows 



£ £ (\N\ + \L\)\5<f>\dvol < AS £ £ (1 ^{ 1+Q rfvol + C £ ^^Ij r + CG[1 + 8a, < 

+ C^j T \l + T)- 1 - 2a E[4 > ](r)dT^ 7 ^£ 2 (l + T)- 7 + 4a g(l + 2a,T)d^j 



Tl 

< AS / - — rr-r rfvol + C / - — \\' dr 
~ J T1 Jvt (l + r) 1+a J T1 (l + r)i+ 2 « 

+ CG[2 + 2a, 1 + 2a] T T \ + CG[l + 2a, a\ T T \ , 
where we used Holder's inequality and Jensen's inequality 

a^W < (1 -7)0 + 76, Va,6>0. 
For the inhomogeneous term |F||<9</>| in (18), we have 

If we choose 

6 °- 100' 

then for all S < So, inequality (18) implies that 

+ G[2 + 2a,l + 2a\Z+G[l + 2a,a]%. 
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Similarly, the energy inequality (9) shows that 



E[a]i r 2 ) < E[a]( Ti > - ■ / j^^-dr + G[2 + 2a,l + 2a\\\ + G[l + 2a, a]^ + D°[F)Z 



We Gronwall's inequality to control the second term on the right hand side of the above inequality. We 
thus have 

E[4>](T2) < E[<f>}( Tl ) + G[2 + 2a,l + 2a]£ + G[l + 2a, a]£ + D a [FY T \. (21) 
Then the above integrated local energy inequality is improved to 

£ L (l + r{i+° dV01 ~ ^ ](Tl) + Da[FVrl + G[2 + ^ 1 + ^ + 0[1 + 2a ' a] "' (22) 

We have shown (21), (22) for all t 2 > t\ > to- We claim that these two inequalities hold for all r 2 > t\ > 0. 
In fact, when n < r 2 < to, the finite speed of propagation for wave equation [21] shows that <f> vanishes 
when r > R = t + R . Hence we can show 



(\N\ + \L\)\d<j>\dxdT < r [ \B^\ 2 dxdr< r E[d>](T)dT, 

Jti Jr<R Jt, 



ti JE T Jti Jr<R 

When considering the energy inequality (9), J T ^ 2?[</>](r)dr can be absorbed by using Gronwall's inequality 
since t 2 < to. Hence we can conclude (21), (22) for all < t\ < r 2 < t - For the case n < t < t 2 , split 
the interval [ti,t 2 ] into [ri,to] and [io7 T 2], 011 which we have two separate inequalities. Combining them 
together, we get (21), (22). Therefore (21), (22) hold for all < n < r 2 . 

We end this section by making a remark. We have used the modified energy flux E[(ft\(r) instead 
of E[4>](t) to make the above argument rigorous. We claim that the inequalities (21), (22) hold if we 
replace E[4>](t) with E[4>](t). In fact, it is sufficient to consider the case when 



%](ti) + D a [F] T T \ + G[2 + 2a, 1 + 2a\ T T \ + G[l + 2a, a] 



is finite. By (21), this shows that E[4>](t) is finite for all t e [ti,t 2 ]. Thus Remark 1 shows that all the 
above statements hold if we replace E[4>](t) with E[4>](t) for r e [ti,t 2 ]. In the sequel, wc no longer use 
the modified energy flux E[4>](t) for the reason argued here. 

3.2 p-weighted Energy inequality 

We revisit the p-weighted energy inequalities developed by M. Dafermos and I. Rodnianski in [5]. Rewrite 
the equation (13) in null coordinates 

-dM + ^ = r{F-N-L), 4>:=rd>, (23) 

where ^ denotes the Laplacian on the sphere with radius r. Multiplying the equation by r p d v ip and 
integration by parts in the region bounded by the two null hypersurfaces 5 Tl , S T2 and the hypersurface 
{r = R}, we obtain 

f r p {d v ip) 2 dvduj + 2 f 2 f rP +1 (F - N - L)d v ipdvdTduj 

JS T , Jti JSr 

r p_1 (p(d v ip) 2 + (2 - p)|yV| 2 ) dvdrdu + [ r p \f^\ 2 dudTduj 

Jill 

= [ r p {d v ^fdvduj + rP (\fyj\ 2 - (d v yj) 2 ) dudr\ r=R . (24) 
Js^, Jt, 



> T 2 

fT 2 
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We claim that we can estimate the boundary terms on r = R as follows 



f 2 r p (IWI 2 - (3„V) 2 ) duodr 

J Tl 



r=R (25) 

< E[<P]{ Tl ) + G[2 + 2a,l + 2a}H + G[l + 2a, a]% + D a [F]? 



JT1 ■ 



Since R is a fixed constant, it suffices to show (25) for p = 0. Thus take p — in the identity (24). 
The energy term on the null hypersurfaces S Tl , S T2 , U£ can be bounded by Efo) + E(t\), which can be 
estimated by using the energy inequality (21). We use the improved integrated local energy estimates for 
to bound the third term in (24). Recall that when r > R, we in fact have the improved lower bound 

1 < /3((l + r)"-l) 1 1 

r~ r{l + r) a (l + r )i+a X 2 J 



instead of ^ 1+ ^i +a we have used in (8) to obtain (18). Thus we actually can show that 

I 2 [ r- 1 \fi)\ 2 dvdudT= [' f ^^-dvol < £7[0](n) + [' f \F - N - L\\B(f>\dvol. 
For the inhomogeneous term, notice that 

f 2 f \r{F - N - L)d v ^\dvdudr < f 2 f \F - N - L\\d<t>\dvo\. 

Jti J S t Jti Js t 

We have already shown that this term can be bounded by the right hand side of (25) in the previous 
section. Thus the inequality (25) follows. 

Now, to make use of the identity (24), we need to control the inhomogeneous term r p+1 (F — N — L)d v ip 
as all the other terms have a positive sign or are bounded. Due to the different structures of F, N, L, we 
discuss them separately. The most difficult term is the linear term r p+1 iV($, <j>)d v ip satisfying the null 
condition. For this term, by Lemma 4, it suffices to estimate the following three terms 

r^ll^VH^!, r p |d$||</>||3W|, r p \d v $\\d^\\d v ^\. 

In application, p <G (0,2). In particular, the coefficients p, 2 — p in (24) are positive. From the decay 
assumptions on $( see Definition 1), we estimate the first term as follows 

2r p |d$||aW>II^Vl < 2r p -5(i + (t - |a;|) + )- i ~ 4c WIM 

C (26) 
<e 1 rP- 1 \d v ^\ 2 + -rP(l + T)- 1 - 8a \d v ^\ 2 1 Vei > 0. 
ei 

The first term will be absorbed if p > ei, 2 — p > a, while the second term will be controlled by using 
Gronwall's inequality. Similarly for the second term r p |9 < i>||0||9 u V'|, we can show 

2r p \d<P\\cj)\\d v iP\ < r p - 1+3a (l + T)- 6a \d v i>\ 2 + r p - 3a (l + t)- 1 - 201 ^] 2 . 

We use interpolation to further bound the first term on the right hand side of the above inequality. Notice 
that 

V ■ > P - 1 + 3a, p < 1 + 2a. 

1 + la 

Using Holder's inequality and Jensen's inequality, we have 

r p - 1+3a {l + T)- ea < (^(I + t)- 1 -") 1 ^ • ((1 + t)-") 1 " 1 ^ <T~ a +r p (l+r)- 1 - a . 
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We use estimate (20) to bound r p 3a <f) 2 . Summarizing, we can show that for p < 1 + 2a 

f 2 f r p \d<I>\\4>\\d v ij\dvdujdT 

1 ^ T2 - (27) 

< rr Q G[0, 0]£ + G[l + a,p] T T \ + G[l + 2a,l + 2< + jf ' JMlL- dr. 

It remains to handle the third term r p |<9^$||<9'0||<9 t ,'i/'|. We estimate this term and the linear term 
r p+1 L(d(f>)d v ilj together due to the similar assumptions on L p (t,x). The difficulty for estimat- 

ing these two terms is that we are not allowed to use Cauchy-Schwartz's inequality as we did previously. 
However notice that the integrated local energy is expected to decay in r( (1 + t)~ 1 ~"). We can put some 
positive weights of r in the integrated local energy such that it is still bounded. To start with, observe 
that when \x\ > R > 1, we have 

r p \d v <$>\\dtp\ +r p+1 \L(d(t>)\ <r p - 1 - 3a \di>\+r p+1 - 1 - 3a \dcj)\ < r p - 3a 



Thus we can bound 

f T2 



f 2 f r p \d v $\\dil)\\d v ip\ +r p+1 \L{d(f>)\\d v xp\dvdujdT < f f r p ~ 3a \d(/)\\d v ip\dvdiudT 

~ h Z Tl ~ a Is (i+t)^ dvo1 + 62 (G[0 ' 0] - } ^ (G[1 + a ^ 2) m 

\ £ rl ~ a L (i+ty+° dvo1 + 62 (G[0 ' 0] - + G[1 + aMl2) 



< 

£2 

for all positive number e2 < 1. Here we have used the fact 

1 — a 

2p - 1 - 5a p < 0, p < 1 + 2a. 

1 + a 

We now have to show that the first term on the right hand side is bounded. We rely on the following 
lemma. 

Lemma 6. Suppose /(t) is smooth. Then for any (3^0, we have the identity 

r s ^f(s)ds=p rv- 1 r f( S )dsdT+T? r 2 /^^. 

J n J ti Jt J n 

Proof. Let 

F(r) - jf* /(*)ds. 

Integration by parts gives the lemma. □ 

Apply the lemma to /? = 1 — a, /(r) = J g (l+^i+c Then the integrated local energy inequality 
(22) implies that 

jf ^ I (l+r{i+° dvo1 ~ £ M( T ) dT + ^"^[^(n) + jf* r-«£> a [F];^r 

+ G[l + 2a, 1 + 2a]£j + G[2a, a]^ + T 1 1_a D a [f]^. 

Since in application only two p-weighted energy inequalities associated to p = 1 and p = 1 + 2a are 
considered, we use interpolation to bound G[2a, a] 

G[2a,a]% < e 2 e 3 G[2a,2a]% + — rf 2a G[0, 0]£ 

£2^3 
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for all positive £3, where 62 is the constant appeared before. 

Our ultimate goal is to derive the decay of the energy flux E[4>](t) on E T . The almost energy flux 
5(0, t) on S T is related to E[<P]{t) by the following lemma. 

Lemma 7. 

E[<P}(r)<g(0,r) + 2 f \d<f>\ 2 + <p 2 dx < E[4>](r). 

Proof. In fact note that 

g(0,r)+2 [ \dcj)\ 2 + 4> 2 dx= [ r 2 (d v $) 2 + d v (r<j) 2 ) + r 2 \f(j)\ 2 dvduj + 2 f \d<j)\ 2 + cj) 2 dx 

Jr<R Js T Jr<R. 

= E[4>}(t)+ / r4> 2 duj + / \d(j)\ 2 + 2<\> 2 dx. 

J lq v T J r<R 

Lemma 1 and Lemma 2 imply that 

5(0, r) + 2 f \dcj > \ 2 + ci> 2 dx<E[4>]{T). 
To prove the other side of the inequality, it suffices to show that 

r<j) 2 {T,R,ijj)dijj < / \d r (/)\ 2 + 2<j) 2 dx. 

r=R Jr<R 

Without loss of generality, assume R>2. Notice that 
r r-R 

?3 



R s [ (f> 2 {T,R,uj)du = [ ( d r (r 3 (j) 2 )dujdr <3 / (f> 2 dx + R [ \d r 4>\ 2 + (f> 2 dx. 

Ju> JO Ju) Jr<R. Jr<R 



Hence 



R [ <f> 2 {T,R,uj)duj < [ \d r (f>\ 2 + 2<f> 2 dx < [ \d<f>\ 2 +2<j) 2 dx. 

Ju Jr<R Jr<R. 

Thus the lemma holds. □ 
Since G[0,0]^ < G[0,0]^, using Lemma 7, we can control G[0,0] in terms of E[4>](t) 

G[0,0]™ <G[0,0]£ < m(r)dr. 

Jri 

Summarizing, we can show that 

r p \d v ^\\dilj\\d v ip\ + r p+1 \L(d<p)\\d v ^\dvdujdT 



< (e 2 + ^ + T A^-\ f 2 E[<j>]( T )dT + rl^Emn) + C T- a D a [F]?dr (28) 
+ T{- a D a [FY T l + -G[l + 2a,l + 2a]% + e 3 G[2a, 2a\\\ + G[l + a,p]%. 

£2 

Here we used the argument in the end of previous section to replace E[<P]{t) with E[(f\(T). We must 
remark here that the implicit constants before the other terms on the right hand side of (28) may also 
depend on e- L . However, since q will be chosen to depend only on R, a, B a P , C , the omitted dependence 
will not affect the argument in the sequel. 
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Finally, we treat the inhomogeneous term r p+1 F ■ d v ip in (24). Since .D 2a [F]£ 2 is expected to decay in 
r, we put some positive weights of r in D 2a [F]^ 2 and estimate it by using Lemma 6 applied to (3 = p— |a, 
/(t) = J St r 1+2a |F| 2 dcr. We can show that 

/ 2 / 2rP +1 F ■ d^dvdrdu 

<e 4 G[p-la,2p-l-2a} T T l + - [ \l + T )P-? a [ \F\ 2 r 1+2a dvol (29) 
2 e 4 J Tl Js r 

< e 4 G[p - |a, 2p - 1 - 2a]™ + — / * tP- 1 "^ 2 "^]; 2 ^ + — r v ~^ a D 2a [F \ T T \ 

2 £4 J Tl £4 

for any < £4 < 1 and r 2 > n > to- 
3.3 Proof of Proposition 1 

Having controlled J^ 2 J g r p+i(p — jy — L)d v ipdvdujdT, we are now able to prove Proposition 1. First let 

Hence for p = 1 or 1 + 2a, the third term in (24) dominates the first term on the right hand side of (26). 
Set p = 1 + 2a in (24) and £4 = 1 in (29). Combining the estimates (25), (26), (27), (28), we infer 
that 

g(l + 2a,T 2 ) + G[0,2a} r T l<g(l + 2a,T 1 )+(e 2 + ^ + ^—] f ' E[<f>]( T )dT + I ' T* a D 2a [F\?dT 

+ Tt- a E[4>}( Tl ) + rl +ha D 2a [FY T \ + £ 3 G[0, 2a]£ + -G[l + \a, 1 + 2a]£. 



£2 2 

Now suppose the implicit constant before £3(5[0, 2a] is Ci, which is independent of £2, £3. Take 

63 = 2^- 

We remark here that we can choose different £j for different values of p. In particular, we conclude that 
esG[0, 2a] can be absorbed by the left hand side. Then apply Gronwall's inequality( Lemma 5). We can 
control the last term j-G[l + |a, 1 + 2a]^ 2 and conclude that 



g(l + 2a,T 2 ) + G[0,2aY T l<g(l + 2a,T 1 ) + T 1 1 + ? a D 2a [F} T T l+ I " 'D 2a [F}? dr 

f J -h a \ / -Ct T--2a\ /-T2 

+ exp f J (^£ 2 + IL_ + J_5_j y E[^(r)dT + ^"^[^(n). 



The integral of the energy on the right hand side can be estimated when we combine (30) with the 
p- weighted energy inequality for p = 1. 

Now take p — 1 in (24). First, we use interpolation to estimate the first term G[l — |a, 1 — 2a] ^ on 
the right hand side of (29) 

G[l - \a, 1 - 2< < (g[1 + \a, 1]%) (G[0,0]^ 2 ) 2a < G[l + 1]™ + G[0,0]%. 

To retrieve the full energy E[4>](t) from <?(0,t), by Lemma 7, add 

r T 2 r 

.12 , ,2, 



in Jr<_R 
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to both sides of (24). Then the integrated local energy estimate (22) restricted to the region r < R and 
Gronwall's inequality imply that 

g(l,r 2 ) + jH E[4>]{T)dr < g(l,n) + (e 2 + ^ + e 4 + ^) jT i?M(r)rfr + £ T- a D 2a [F]?dr 

+ ri-^^Kn) + rl- a D 2a [F]H + rf 2tt G[0, 2a]£ + -G[l + 2a,l + 2a]£, 

where we choose £3 = 1. Assume the implicit constant before J^ 2 E[<j>](T)d,T in the above inequality is 
C3, which is independent of e 2 and €4. Then take 



64 -2C 3 - 



We get 



9(1, r 2 ) + jH £[<A](r)dr < 5 (1, n) + (e 2 + ^ + ^) jf" S[0](r)dT + jH T- a D 2a [F]?dr 

+ Tl- a EmT 1 ) + Tl- a D 2a [F} T / i +T^ a G[0 7 2a}ll + -G[l + 2a^ + 2a}^ i . (31) 

£2 

Now let C 4 be the implicit constant before J^ 2 E[<f>](T)dT in both (30) and (31), which is independent 
of e 2 . Then let 

1 

62 ~acI 



and choose a constant To > io such that 
In particular, for t\ > To, we have 



T-J« < " 



/ — a — 2a \ / rp — a rp — 2a\ -i 

We combine (30) and (31) together to control J^ 2 E[<P\(t)(It. For r 2 > n > T , we first estimate 
G[0, 2a]^, G[l + 2a, 1 + 2a]£[ in (31) by using (30). Then combining all them together, we can show 
that the coefficient of J^ 2 E[<j)](T)dT on the right hand side can be bounded by 

1 „ ^-2a e C 4 T- 2a e 3 
2 +C ^ 2 + ^^2<4- 

Thus J^ 2 S[(^](r)dr can be absorbed and we can conclude that 

5(1,7-2) + P BM(r)dr < 5(l,n) + T{~ a D 2a [F]^ + E[<t>]( Tl ) + P r- a D 2a [F\?dr 
J n «/ti 

+ rf 2a H r^C 2 «[F]; 2 dr + rf 2Q . 9 (l + 2a, n), 

which, in turn, improves (30) to 

g(l + 2a,r 2 ) + G[0,2 a }ll<g(l + 2a,T 1 ) + rl- a EmT 1 )+rl + ^ a D 2a [FY T l + P r^D 2a {F}?dT. 

Jti 

This proves (16) and (17) for all t 2 > r x > T . 
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For t < T\ < t 2 < T , we make use of the boundedness of r. Let e 2 = 1. Inequality (30) shows that 
g(l + 2a,T 2 ) + G[0,2a} r T l<g(l + 2a,T 1 )+ I ' * S[0](r)dr + £> 2tt [F]£ + E^n). 

Jtl 

By Lemma 7, we have 

G[l + 2a, a]% < G[0, 2< + G[0, 0]£ < g(l + 2a, n ) + P E^dr + D 2 "[F]£ + B[^](ti). 

Jtl 

Combining with the energy inequality (21), we obtain 

EM-m) S 5(1 + 2a, n) + f * £[0](r)dT + L> 2a [i^ + £7[0](n). 
Thus Gronwall's inequality indicates that 

£[0](r)dr < 5 (1 + 2a, n) + D 2a [F\% + E[cb]( Tl ) 



I 

Jtl 



as Ti < r 2 < T . Hence (16) and (17) follow from (30), (31). 

For n < T2 < £o> the finite speed of propagation for wave equation [21] shows that g{p,r) vanishes. 
Thus (16), (17) hold. For general r 2 > T\ > 0, divide the interval [ti,t 2 ] into three(possibly two) 
such intervals: [ri,t ], [io,^o] an d [T ,t 2 ]. Then (16), (17) follow by combining those three(or two) 
inequalities together. This completes the proof for (16), (17). 

Having proven (16) and (17), we can improve the integrated local energy inequality (21) and the 
energy inequality (22) as follows: Integrate (17) from n to r 2 . We obtain 

G[2 + 2a, 1 + 2a) T T \ + G[l + 2a, a\ r T \ < G[2 + 2a, 1 + 2a] T T \ + (1 + n)" 1 "^^, 2a] T T \ 
< E[<f>](n) + r^ a D 2a [F] T T 2 i + (1 + n)- 1 - 2 " (g(l + 2a, r 2 ) + £ T^D 2a {F}?d^J , 

which, together with (21), (22), implies (14), (15). We thus finished proof for Proposition 1. 
To show Corollary 2, take p = 1 + 2a in (20). Interpolation shows that 

r 2a 4> 2 dvdu <E[4>]{t) +E[4>]{T) 1 -^g(l + 2a, T )^. 



L 



Using Jensen's inequality, we have 

rT 2 r „2ai2 



Therefore for r 2 > n >t , Lemma 4 and Proposition 1 imply that 
D 2a [N] T T \ = ( 2 f \B af) d a § ■ dp<f)\ 2 {l + r) 1+2a dvol 



(l +r )l+a (! + T )l+8a (l + r )l+4a (1 + r )l+8a 



ri 

~ X. (TT7F^ V01 + ^ + **' 2a] " + /J (l+ [ r) ( i+^ dT + ° [2 + °' 1 + 2a] 

< EMn) + D«[FY T \ + (1 + Tl )- a D 2a [F}ll 

+ (l + r 1 )- 1 - 2a (g(l + 2a,T 2 )+ J\i a D 2a lF}?dT^j . 
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For ti < t 2 < t a , notice that 

D 2a [N]H< f 2 E[4>}{r)dr. 
For the linear terms L(d<j)), we can show 



D 2a [L] T T \ = r f \L(dc/>)\ 2 (l + r) 1+2a dvo\< f (1 + r y 2 - 6a \d<f>\ 2 (l + r) 1+2a dvol 



< i \d(j)\ 2 

(1 + r) 

The corollary then follows from (14) and (15) 



^ / / 7Z — rfvol. 



4 Decay of the Solution 

Under appropriate assumptions on the inhomogeneous term F, Proposition 1 leads to the decay of the 
energy flux E[<p](r). After commuting the equation with the vector fields Z, we obtain the pointwise 
decay of the solution outside the cylinder {(t, x)\\x\ < R} by using Sobolev embedding and inside the 
cylinder by using elliptic estimates. 

Proposition 2. Suppose there is a constant C\ such that 

D 2a [FY T \<C l {l + r 1 )- l - a , VT 2 >n>0. 

Then for solution <p of the linear wave equation (13), we have energy flux decay 

^(rJ^^oTCOa+r)- 1 - 

Proof. Since the initial data are supported in the region {|x| < i? < R}, the finite speed of propagation 
shows that g(l + 2a, 0) vanishes. Take t\ = in (17). We get 

g(l + 2a,T)= [ r 1+2a (d v iP) 2 dvdLo<Ci+e 2 E a (32) 



and 



/ 2 / r 2a (d v ip) 2 dvdudT < G[0, 2a] T Q 2 <d+ e 2 E . (33) 



In JS T 

We claim that we can choose a dyadic sequence {r„ — > oo} such that 

/ r 2a (d v ^) 2 dvdw < (1 + t^- 1 (d + e 2 E ) , (34) 

where r„ satisfies the inequality 7~ 2 r„ < r n _i < 7 2 t„ for some large constant 7. In fact, there exists 
r n € [7",7 n+1 ] such that (34) holds. Otherwise 



/ / r 2a (9 t ,V) 2 ^wrfr>ln7(e 2 J B + Ci), 



which contradicts to (33) if 7 is large enough. 

Take r = r„ in (32). Interpolate with (34). We obtain 



/ r(c^) 2 cfcdw < (1 + r n )- 2a (e 2 E Q + &) . 
Js T „ 
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Then the inequality (16) implies that for r > r„ 

f r(d v ^) 2 dvdLo + f E[4>](s)ds < (1 + r„)- 2 « (e 2 E + d) + ri- Q i?[^](T n ). (35) 

On the other hand the energy inequality (15) shows that for all s < r 

M < ^M(s) + (1 + a)" 1 " (e 2 E + Ci) . 

In particular 

BM(ti) < £?M(0) + e 2 £ + d < e 2 £ + d- 

By (35), we have 

(r - r n )£W](T) - J\\ + .s)- 1 -" (e 2 £ + d) ds < (1 + r„)- 2 « ( £ 2 £ + d) + r^"^]^). (36) 
In particular for n = 1 

Let r = r„+i in (36). We obtain 

(t„+i - r»)E[ir„ +1 ) < (1 + r„)" Q (e 2 £ + d) . 
Since r„ are dyadic, we have 

E[<j ) ](T n )<T- 1 - a (e 2 Eo + C 1 ), Vn. 
Finally, for t g [r„,r„ + i], we can show 

%](t) < ^M(r„) + (1 + Tn)- 1 -" (e 2 E + d) < (1 + r„)- 1 " a (e 2 £ + d) < (1 + r)" 1 "" (e 2 £ + d) . 

□ 

With the energy flux decay, we can obtain the decay of the spherical average of the solution. 
Corollary 3. Assume that there is a constant C\ such that 

D 2a [F] T T l KCiil + n)- 1 -", Vr 2 >T!>0. 
Then on the hyper surface S T , we have 



[ \r<j>\ 2 diu<e 2 E Q + C u r > R, 

[ r\^\ 2 du: < (1 + r)- 1 -" (e 2 E + d) , r > R. 



Proof. By Proposition 2, the first inequality follows from (19) and (32). The second one follows from 
Lemma 1. □ 

In order to obtain the pointwise decay of the solution which is usually a consequence of Sobolev 
embedding, we need energy estimates for the derivative of the solution. For this purpose, we commute 
the equation with the vector fields and T. Under appropriate assumptions on the inhomogeneous term 
F, we hope to derive the same energy decay for VL k T^<p. Denote 

N{<h.,<h) = B^d^-dvfo, V0i,02 G C°° (M 3+1 ), 
where we recall that the constants B a @ satisfy the null condition. 
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Lemma 8. Let Z be SI or T. Then 



Proof. Notice that [SI, d r ] = [SI, d t ] — [SI, )t7] = 0. The lemma then follows from the fact that B a P satisfy 
the null condition. □ 

Based on Corollary 2, we are able to prove the decay of the energy flux of Z^cj) after commuting the 
linear equation (13) with Z^ . 

Proposition 3. Assume that there is a constant C\ such that the inhomogeneous term F in (13) satisfies 
the following condition 

D 2a [Z tj FY T \ < d(l + n)- 1 -", Vr 2 > n > 0, V/3 < [3 

for some multiple indices \(3 \ < 4. Assume $, L^(t,x) satisfy the conditions in Theorem 2. Then we 
have 

E[ZU]{r) < (Ci + e 2 E ) (1 + n)" 1 "", (37) 

D 2a [N(Z^<S>, Z^Wrl + D 2a [Z^L» ■ Z%4>)]1^ < (d + e 2 E ) (1 + n)" 1 " (38) 

forVf3</3 , <4. 

Proof. We prove the proposition by induction. When (3 — 0, (37) follows from Proposition 2. Since 
Z^ 1 ^ is (S, a, to, R\, Co)-weak wave, V|/?i| < 4, Corollary 2 and inequality (32) imply that 

D 2a [N(Z^<S>, <j>)Y T \ + D 2a [Z^L» ■ 8^]% < (Ci + e 2 E )(l + n)" 1 "", V|ft| < 4. 

Assume that (37), (38) hold for all (3' < (3. Commute the equation (13) with Z$ . Using Lemma 8, we 
have the equation for Z^<j) 

U{Z fi 4,) + N{^,Z f} cf))+L{Z fi cj))^Z l3 F~ N{Z h §,Zl )2 <f>) + Z^L" ■ Z h d^. (39) 

/3i+/32</3,/3 2 </3 



Since @2 < ft, by the induction assumptions, we get 

Z f3 F - N{Z f3l $, Z l)2 cj)) + Z^U* ■ Z^d^ 

02<P 



< (C 1 +e 2 E )(l + T 1 ) 



-l-a 



Hence for Z^<j>, inequality (37) follows from Proposition 2 and inequality (38) follows from Corollary 2 
and Proposition 1. □ 



Since the angular momentum SI is vanishing for r = 0, we are not able to obtain the pointwise bound 
of the solution in the cylinder {|a;| < R} by commuting the equation with SI. We instead rely on elliptic 
estimates and the vector T — d t as commutators. 

Lemma 9. Assume that there is a constant C\ such that 

D 2a [F] T T \ + D 2a [d t FY T \ <C 1 (1 + ri)- 1 -", Vr 2 > n > 0. 
Then for solution of the linear wave equation (13) ? we have 

3 

f \d 2 (t>\ 2 dx^ V f Id^dx^tEoe' + C^il + T)- 1 -*. 
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Proof. We first assume that r > R. Take /3 = (0, 1) in Proposition 3. We have 

E[T^](n) + £ £ { l d + r f +a dxdT + D 2a [T^N(^ <P) + T'L{d<f>)]^ < (E e 2 + &) (1 + n)" 1 " 
for all j < 1. Using elliptic estimates [6] , we can show that 

/ \d 2 <f>\ 2 dx = ^2 \d tJ cl)\ 2 dx + 2^2 \d a d t <j)\ 2 dx 

< E[d t 4>](r) + f |A0| 2 + \4>\ 2 dx 

Jr<2R 

<E[d t 4>]{ T )+ ( \d tt <j> + F-N{<S>,<j>)-L{d<f>)\ 2 + <j> 2 dx (40) 

Jr<2R 

< E[d t <f>](r) + / \d tt <t>\ 2 dx + V / / |7*>| 2 + |T>| 2 + \T> N\ 2 + \T>L\ 2 dxdt. 

Jr<2R j<l^ T Jr<2R 



Consider the region bounded by X t -.r and t — t. Take X = T in (7). Note that the vector field T is 
killing, that is K T [d t cj)] — 0. We can conclude that 

[ jJ^Kda = f jJ[d^Kda+ f f {TN + TL-d t F)d tt <j>dvol. 

Jr<2R JE T - H n{i<T} JT-RJr<R+t-T 

Apply Cauchy-Schwartz inequality to the last term. We obtain 
/ I \d t F -TN -TL\\d tt <j)\dvol 

Jt-R Jr<R+t-T 

< [ \dtF -TN- TL\\l + rr+* + dxdt 

-I f n'f'tL dxdt + D2a & F Vr- R + D 2a [TNY T _ R + J 2 "[T£];_ fl 

Jr-RJ^t U + r J 

Hence we can estimate 

/ \d tt ^\ 2 dx<2[ jT[d t 4>}n»da<E[d t (f>}(T-R)+ f f \d t F - TN - T L\\d tt <f>\dvol 

Jr<2R Jr<2R JT-RJr<R+t-T 

< {E e 2 + C 1 ) (l + n)- 1 -". 

Then from (40), we get 

/' \d 2 4>\ 2 dx < (E e 2 + d) (l + n)" 1 -" + ]T f T+R I n f 1i+a + l TJF l 2 + l TJiV l 2 + l T ^l 2 

< (s e 2 + d) (i + n)- 1 -" + ^£> Q [r j 'F];tg + i? 2Q [T j iv];;^ + i? 2a [r j L];;g 

< (E e 2 + C 1 ) (l + n)- 1 -". 

Thus we have proven the lemma for t > R. When t < R, the finite speed of propagation shows that the 
solution of (13) vanishes when \x\ > t + Rq. Thus wc can replace r — R with in the above argument. 
And the lemma still holds. □ 

A corollary of the above lemma is the following pointwise decay of the solution in the cylinder {r < R}. 
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Corollary 4. Assume that there is a constant C\ such that 

D 2a [F] T T l + D 2a [d t F]H < Ci(l + n)- 1 -", Vr 2 > n > 0. 
TTierc /or solution <p of (13), we /icwe 

^^(d+e^a + r)- 1 -", r<i?. 
Proof. Using Sobolev embedding and Lemma 9, when \x\ < R, we can estimate 

3 

^ 2 < / y i5„0i 2 + ^ 



'r<i? Jt Jr<R 



<{C 1 +e 2 E ) (1 + t)- 

where the last step follows from the integrated local energy inequality (14) restricted to the region 
r<R. □ 

5 Boostrap Argument 

To solve our nonlinear problem (4), we use the standard Picard iteration process. We prove, by a 
boostrap argument, that the nonlinear term D 2a [F]T r l decays, which leads to the decay of the solution </>. 
We still denote the quadratic nonlinearity A afi d a <pdp4> of F in (4) as N{<p,<p) — A^d^d^cj), in which 
the constants A^ v satisfy the null condition. 

Proposition 4. Suppose Z k $> is (5, a, to, Ri,Ci)—weak wave for all \k\ < 3. Assume 

\d 2 Z^\<C u V|/3|<1. 

Assume the functions L^(t,x), h^it.x) satisfy the conditions in Theorem 2. If the nonlinearity F in 
(4) satisfies 

D 2a [Z?F]H < 2E e 2 (l + n)" 1 "", V|/3| < 3, Vr 2 > n > 0, 
f \\7Z !3 F\ 2 dx < 2E e 2 (l + T)- 1 - a , V|/3| < 1, Vt > 0, 

Jr< 



lr<R 



then 



D ^ [ ZPFY T l<E 2 e i {l + T 1 )- 1 - a , V|/3|<3, Vr 2 > n > 0, (41) 
/ iVZ^Ffdx^Eyil+r)- 1 - , V|/3|<1, Vr>0. (42) 

Jr<R. 

Remark 9. If the given function $ is assumed as in Theorem 2, that is, Z 13 ® is (S 7 a 7 t 07 Ri,Ci) — weak 
wave for all \(3\ < 4 and \d 2 Z^^\ < C\ for all |/3| < 2, then the above proposition holds if we replace 
\P\ < 3, |/3| < 1 with \j3\ < 4 7 \j3\ < 2 respectively. The reason that we formulate the proposition as above 
is that three derivatives are the minimum to close the boostrap argument. Four derivatives is needed to 
obtain C 2 solution of the equation (4). 

The proof for Proposition 4 is quite similar to that in [22]. For completeness, we repeat it here. Since 
higher order nonlinearity decays much better, we only consider the quadratic nonlinearities N{(j>, </>) and 
Q(d(j)) = h tiV {t 1 x)d^4>d v (j). First, Lemma 8 and the assumptions on h^ v (t,x) 

\Z !i h^\ < (l + r)-i a , V|^|<4 
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imply that 

D 2a [Z f) FY T \< [ 2 I \N{Z fll (j>,Z 02 <j>)\ 2 (l + r) 1+2a + \dZP 1 <j>\ 2 \dZfc\ 2 (l + r) 1 - a dxdT 

< r f \ 9 <h\ 2 \ d <h\ 2 dxdT + ]T r f \ N (^^)\ 2 r 1+2a dxdr (43) 

01+P2<P Jt1 Jr - R Pl+P2<0 Tl St 

+ E I \d<h\ 2 \d<h\ 2 {^ + rf- a dxdT, 

Pi+I3 2 <P Jti JS ^ 

where we denote fa = Z^ 1 ^, fa — Z^ 2 <j). We estimate the three integrals on the right hand side of (43) 
separately. We use elliptic estimates as well as the extra boostrap assumption (41) to bound the first 
integral. Estimates of the second integral rely on the null structure of N(fa,fa). The third integral 
follows from the integrated local energy inequality (14). 

5.1 Proof of the (41) in the region r < R 

When r < R, we use elliptic estimates to obtain the pointwise bound of the solution. However since one 
can only obtain elliptic estimates in a smaller region, we divide the region r < R into two parts: r < -j 
and r > -j. In the smaller region r < § , we use elliptic estimates while in the region r > |, we use 
Sobolev embedding. 

Recall that fa = Z^fa fa = Z^fa \fa\ + \/3 2 \ < 3. Without loss of generality, assume \ fa\ < \/3 2 \. In 
particular we have \fa \ < 1. For r < we claim that 

\dfa\ 2 <E e 2 (l + r)- 1 - a , r<|. (44) 

We first verify (44) for dtfa- By (39), dtfa satisfies the following equation 

U(d t fa)+N{<S>,d t fa) + L{d t fa)=F l . 
Since \ fa + (0, 2)| < 3, estimates (38) imply that 

£> 2a [iq]^ + D^F^ < E e 2 (l + n)" 1 "". 
Thus by Corollary 4, we have 

\d t fa\ 2 <E e 2 (l + r)- 1 - a , r<|. 

For Vfa, notice that fa = Z^fa \fa\ < 1 and \d 2 Z^<S>\ < 1, \dZ^L^\ < 1. Using elliptic estimates and 
Lemma 9, we have for \x\ < ^ 

HV4>i||V < / fl \d^fa\ 2 + \Vfa\ 2 dx 

< f IVA^P + IV^I 2 ^ 

Jr<R 

< I |V (dufa + Z^F - Z^N($,fa - Z^L(dfa) | 2 + \Vfa\ 2 dx 

Jr<R 

< E[d tt fa](r) + E[fa](r) + f \VZ^F\ 2 + \V 2 fa\ 2 + |V 2 0| 2 + \Vfa\ 2 + \V0\ 2 dx 

Jr<R 

<E e 2 (l + T 1 )- 1 - a , 

(45) 
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where Bi R denotes the ball {r < \R} in M 3 . Hence we have proven (44), which implies that 
/ 2 / |d0i| 2 |c^ 2 | 2 cfa;< ( 2 {l + ry^Eoe 2 f \dfo\ 2 dxdr 

J Tl Jr<§ Vti Jr<§ 

<E e 2 r{l+r)- l - a E[^]{T)dT 



IT1 

^2 4/-, , _\-l-a 



In the region -j < r < R, we use the angular momentum SI. By Sobolev embedding on the unit 
sphere, we have 

f \d<j>i\ 2 ■ \d<f> 2 \ 2 dco < J2 I \d<h.'?du- f \d(t> 2 \ 2 duj, (46) 

J UJ y J UJ J UJ 

where we still denote <p v = Z^'<t> for /3 V < fix + (2, 0). Notice that \/3 v \ + \/3 2 \ < 3 + 2 = 5. Without loss 
of generality, we assume | 1 < 2. Thus by Lemma 9, we have 



[\d<p 2 \ 2 du< f \d<p 2 \ 2 + \d r d<p 2 \ 2 dx<E[4> 2 ](T)+ f ^^dx^ii + ry^Eo 

Jul J Jr<R 



^, (47) 

lr<R 

l<r<R 



where we have used (38), (39) and the assumption \ fi 2 \ < 2 to verify the conditions in Lemma 9. Since 
\Pv | < 3, we can show that 



/ 2 / {dfafldfol 2 dxdr < [ 2 [ f Idfafld&l 2 duj r 2 drdr 

Jt! J±R<r<R Jn J§ Jul 

-zZ[ 2 [ [ \ d <f>i'\ 2 du [ \dcj) 2 \ 2 duj r 2 drdr 
< V rV + r)- 1 -"^ 2 / \d<f> v \ 2 dxd T 

v Jn Jr<R 



<E 2 e\l + r)- l - a . 

Summarizing, we have shown 



V r [ Id^ld^dxdT^E^il+T)-^ 



/3i+/32</3 



Remark 10. We remark here that (47) is only true when r is bigger that a constant. That is why we 
need to distinguish the two cases r < \R and r > \R. 

5.2 Proof of (42) 

Note that when |/3| < 1 and /3i + fi 2 = (3, we have (3\ — or (3 2 = 0. By (46), we have 
f \VZ p F\ 2 du< f \dZ^(f>\ 2 -\d 2 <f>\ 2 + \d(f>\ 2 \d 2 Z^\ 2 duj 

J UJ J UJ 

~ /2 I \9Z^\ 2 duj- f \d 2 Z fi 'c^\ 2 duo+ f \dZ p '(j>\ 2 du- f \d 2 Z P (j)\ 2 du, 

\j3'\< 2 

where as pointed out previously, we only have to consider the quadratic nonlinearities N((p,(p), Q{d(f). 
For r < i.R, the inequality (44) shows that 

\dZ^\ <e 2 E (l + r)- 1 - a , V|/3|<1. 
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For |i? <r<R, the inequality (47) implies that 

f \dZP'cj>\ 2 ( Lj<e 2 Eo(l + T)- 1 - a , V|/3'|<2. 
On the other hand, using Lemma 9, we obtain 

/ \&Z< i '<l>\ 2 dx<e i E (l + T)- 1 - a , V|/3'|<2. 

Jr<R 

Therefore, for all \ j3\ < 1, we can estimate 

/ \VZ^F\ 2 dx<e 2 E (l + r)- 1 - a V / \d 2 Z^' ( t>\ 2 dx<Ele\l + T)- 1 - a . 

Jr<R \f3'\<2 Jr ^ R 

Hence we have proven (42). 

5.3 Proof of (41) in the region r > R 

We first consider the quadratic term N(<f>i, (fo). The p- weighted energy inequality is about tp — rcj) instead 
of <p. For this reason, we expand N(<j>\,<j>2) in terms of i[>. 

Lemma 10. Suppose 7V(0i,0 2 ) = A a P d a <t>idp<l)2 with constants A a/3 satisfying the null condition. Then 
r 4 \N(<j> u <j> 2 )\ 2 < fy 2 + <j> 2 ■ r 2 d 2 r <j, 2 + |yVi| 2 |W 2 | 2 + \d v H 2 \duH 2 + |^i| 2 |S^| 2 , (48) 
where v = , u = . 
Proof. In fact, notice that 

r 2 N((f> U (f>2) = <pi<h+r(<i>i<h)r + N(ipi,ip2) 

and 

MihM < \d v H ■ \d u H + \d u tpi\ ■ \d v ip 2 \ + \fH ■ \fH- 

Hence the lemma holds. □ 

To estimate the second term in (43), it suffices to handle the terms on the right hand side of (48). 
We estimate the first three terms in a uniform way. Let $i be <pi or ^Vii ^2 be 2 , rd r <p2 and ^7tp2 
respectively. Recall that <f>\ = Z 131 ^, </> 2 = Z 132 ^, + \f3 2 \ < 3. Using Sobolev embedding on the unit 
sphere, we have 



/ \$i\ 2 \$2\ 2 dLu<Y^ I \$v\ 2 duj- f |$ 2 '| 2 ^, 

J U) y 2 ! 



where we let 

Uv < ft + (2,0), fc=fa, if|A|<l, f4Q) 
\Pv=Px, < ft + (2,0), if|ft|<l. 

In particular |ft<| + |ft>| < 5. For the third case when $i = y^ij $2 = ^i>2, without loss of generality, 
we assume |ft/| < 2. Since y^i = ^01, $i' can always be written as Z^4> for some \fi\ < 3. Thus by 
Corollary 3, we have 

r 2 / \^ v \ 2 duj < e 2 E , r>R. 

J UJ 

Recall that $ 2 ' = <f>2 f , rd r <j)2' or ^^2'- We always have 



\bm 2 

(i + r) 3 +« ~ (i + r y+ a ' 



< 
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Then the integrated energy inequality (14) implies that 

PT2 P PT2 POO P 

/ / r 2Q - 3 $?$2 dvol = III r^-^lQldvdudT 

pT2 POO P P 

~S / / r 2a - 3 r 2 / l^i'l 2 ^ / |* 2 '| 2 dwdt;dT 

p 2 , Jn JS T r 

< e 2 soV f 2 / ^ l<&2 ;' 2 + rfvoi 

where we recall that a < We hence have estimated the first three terms in (48). 

It remains to handle the last two terms \d v ipi\ 2 \d u ip2\ 2 , \d u ipi\ 2 \dvip2\ 2 ■ Since they are symmetric, it 
suffices to consider \d v ipi\ 2 \d u ip 2 \ 2 . Recall that Vi = rZ t3l (f>, ip 2 = rZ fj2 (j), + |/3 2 | < 3. Define /3i',ft>' 
as in (49). In particular \f3y \ + \02' | < 5. We have two cases according to ft, i — 1, 2. 

We first consider the case when |ft'| < 2. The idea is that we bound uniformly and then 

control IftiV^I 2 by the energy flux through the null hypersurface v — constant. The following lemma 
shows that the energy flux through v — constant is bounded. 

Lemma 11. Consider the region D — [m, u 2 ] x [v\, oo) cS r x [n, t 2 ]. Under the conditions of proposition 
4, we have the energy flux estimate through the hypersurface v — const 

I 2 / ' {d u ^ 2 ) 2 dudio<e 2 E (l + T 1 )- 1 - a 1 

where ip 2 = r<p 2 = rZ^cj). 

Proof. Back to the energy equation (7), take X = T on the region D. We have 

J 2 jT[(f> 2 }n»do- + JJ jT[<j> 2 ]ni*do- = J J Jl^n^da + J J^^n^da 

tf>^l,U-K! V>V-i,U—U2 1 



+ / U(p 2 ■ dtfodvol. 

ID 



Using the estimates (38), we conclude that 



D 2a p<h]% <E e 2 (l + r)- 1 - a . 
Then by the integrated local energy inequality (14) and the energy inequality (15), we can show that 

/ / r 2 (d u fo) 2 du;du < 2 / J^^da < E e 2 (l + n)" 1 "", 

J Ul J Ul J Ul 

where notice that D C S T x [ti,t 2 ]. Thus by Corollary 3, we get 

<E e 2 (l + r 1 )- 1 - a . 



U2 f rU2 r r 

{d u ip2) 2 diodu = / / r 2 (d u (j>2) 2 dudu + / rcjy^dio 

J Ul J Ul J Ul 



ui 

□ 
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We continue our proof of (41) when \(3i> \ < 2. Lemma 11 and Sobolev embedding on the unit sphere 
iw that. 



show that 



In Js 



CO fu{v) 



1 / \d u ip2\ 2 \d v ipi\ 2 dujdudv 

J v T1 J U T ^ J UJ 

poo pu(v) [• 

< / \d u ip2'\ 2 dujdu • sup / r 2a ~ 1 \d v / ipr\ 2 dcu dv 
<e 2 S (l + T 1 )- 1 - a f sup f r^-^dyipvfdujdv 

Jv T1 u Juj 

< e 2 E (l + n)- 1 - a f sup f r 2a \d v iJv\ 2 d^dv, 

Jv T1 u Juj 

where Pi> , /3 2 > arc defined as in (49). Now, for all u <G [u T1 ,u(v)], we have 
r 2a (d v ^v) 2 < r 2a (d v M 2 \ u=ui + / r 2a {d v ^ v ) 2 du 



r 2a {d u d v ^ v fdu+ / r^-Hdv^vfdu 

J Ui 

< r 2a (d v tp v ) 2 \ u=ui + / r 2a (dM 2 du 

+ / r 2a (^ v ) 2 + r 2a+2 \Z^' {F - N - L)\ 2 du, 

where we use the wave equation (23) in the last step and u\ — u T ,u 2 — u(v). Integrate on the unit 
sphere. We obtain 

f sup f r 2a \d v ipy\ 2 dcudv < f r 2a {d v ip v ) 2 dvdu 

+ f 2 i r 2a {d v ibv) 2 + (fnfa,) 2 r 2a +r 2a+2 \Z^' (F - N - L)\ 2 dvdLodr, 
where note that f = 

We claim that the above inequality can be bounded by e 2 _E (up to a constant). In fact, the first term 
can be bounded by e 2 E by (32); the second term can be bounded by e 2 E by (33); the third term can 
be controlled by (1 + Ti)~ 1 ~ a e 2 E by the integrated local energy inequality (14)(notice that \f3y \ + 1 < 3 
and 2a < 1 — a as en < |); the last term can be estimated as 



-l-a 



D 2a [Z^' F] T T \ + D 2a [Z^' N] T T \ + D 2a \Z^' L\ T T \ < E e 2 (l + n ) 
where wc use the inequality (38). Summarizing, we have shown 

f sup f r 2a \d v tp 1 \ 2 dudv<e 2 E . 

Jv T1 u Ju, 

In particular, for fixed r > R, we have 

[ 1 f \d v ^\ 2 (t,r,u)dLjdt < e 2 E , i/>i = rZ^, \/3\ < 2. (50) 

Therefore 



r r2a ~ 3 I |9^i| 2 |9„^ 2 | 2 dvoi < E 2 e\\ + Tl y 



■l-a 
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When \p2'\ < 2, we control Ic^il 2 by the energy and bound d u ip2 uniformly Similarly, using (32), 
(33) and Sovolev embedding, we obtain 



1-1 JS T 
T2 



r^-^dyip^ldu^dvdujdT 



< f* f r^-^d^jyl 2 ■ f \d u il) 2 >\ 2 dLO dvdujdr 

< [ 2 [ r^-^d^yl 2 ■ r - a f \d u tp 2 >\ 2 duj dvdudr 

J T\ J St J OJ 

<e 2 E [ supr- a f \d u ^2>\ 2 dudT, 

Jti v Ju 

where 0i',02' ar e defined as in (49). For all v, we have 



(51) 



r- a (d u ^f<r- a (d u ^) 2 \ 



+ 2 



/ r~ a \d u i 

J V 



r-^du^dv 
ip2< ■ d v d u tp 2 '\dv 
< r- a {d u ^f\ v=VT + / r- 1 - a \8M 2 dv 

T2 JVt 
pOG />OG 

+ I r- 1 - a (d u ^) 2 dv+ I r l - a {d v d u ih) 2 dv 

dv 



< r- a {d^2>)\ =VT2 + 



{d u ^ 2 



-1+a 



+ / r 1 - a (^ 2 ') 2 dv+ / r i -°-\Z^'{F - N - L)\ 2 dv 7 
where v T2 = R ~\ 2 T2 . Integrate on the unit sphere. We get 

[ T2 supr- a [ \d u ^ 2 '\ 2 dudr< H [ r~ a (d u ^) 2 \ v=VT dr 

J Tl V J Ul J Tl J U) T2 

\2 j_ J-o^zP* (F-N- L)\ 2 dvdiudT. 



+ 



We claim that it can be bounded by (1 + n) 1 a e 2 E up to a constant. In fact, the first term can be 
estimated as follows 

I 2 I r- a (d u ^2>) 2 \ v=VT dr= f ^ f r- a \dM 2 dudu< f " f \d u ^ 2 > \ 2 dujdu < (1 + T l y 1 - a e 2 E a 

J Ti J U T2 J U T1 Ju J U T1 Ju 

by Lemma 11; the second and third term can also be controlled by (1 + T 1 )~ 1 ~ a e 2 E by the integrated 
local energy estimates (14)(notice that |/?2'| < 2); the last term can be estimated as follows 

D 2a [Z^'F] T T \ + D 2a [Z^' N] T T l +D 2a [Z^'L] T T \ < E e 2 (l + n)" 1 -" 

by the inequality (38). Hence 

f 2 supr- a f \d u ^2>\ 2 dujdT<Eoe 2 (l + T 1 y 1 - a . (52) 

Jti v Ju 

Plug this into (51). We obtain 

P f r 2a - 3 \d v ^\ 2 \d u ^2\ 2 dyol < e 4 E 2 (l + n)- 1 ^. 

JT! JS r 



29 



Therefore using lemma 10, we have shown that 

[ T2 f |7V(0 1 ,0 2 )| 2 r 1 + 2 «dvol<e 4 J B o 2 (l + n)- 1 - a . 

To show (41), it remains to estimate the quadratic term Z l3 Q{d<j)). By (43), it suffices to consider the 
third integral in (43). Notice that 

\d<h\\d<h\ < r- 2 \d v MduH + £ \Z^'4>\\dZ^cl>\ + (\Z<f>\ + \d v <t>\)\az 3 <t>\. 

\Pl>\ + \P2>\<4 
Pl>,02><3 

The first term has already been estimated considering that r 1_a < r 1+2a (on S T , r > R > 1). For the 
second term, if \(iv\ < 1 or \f3 2 >\ < 1, then using Sobolev inequality on the unit sphere, we have 

! \M 2 \dM 2 <^< I \<t>r>\ 2 cLj ■ / \d<i>2»\ 2 dLo < r- 2 e 2 E £ f \d<j> 2 „\ 2 (ko 

J " |/3i»|,|/32''l<3 \P 2 »\<3 Ju 

by (19) and (32). Then using the integrated local energy inequality (14), we obtain 

fj j s \d<p 1 \ 2 \dM 2 (l + r) 1 - a dvol<e 2 E £ 2 J s dvo\ < e*E 2 (l + n )~^. 

If both | p v | > 2, \p 2 , | > 2, recall that \[3 V \ + \/3 2 >\ < + 1 < 4. We conclude that \/3 v \ = \/3 2 >\ = 2, for 
which we use the embedding 



I \Z 2 ^\ 2 \dZ 2 4>\ 2 d^ < \\Z 2 <pf Hl(s2) \\dZ 2 4>\\% 1{s2) < r- 2 e 2 E \\dZ 2 cb\\ 2 Hl 
I" [ \d^\ 2 \d^\ 2 {l + rf-d,o\< Y^Eo I" I l^gg rfvol^e^a + n)- 1 -" 

Jti J S T \8\<1 1 T 



l (s 2 )- 

Hence 

' ' ^ e 2 E f \ dZ2Qf! ^ 

For the third term (\Z<j>\ + \d v (f)\)\dZ 3 (f>\, using the integrated local energy inequality, it suffices to show 
that 

(|Z0| + |cU>|) 2 (l+r) 2 <e 2 E Q . 
In fact, by Corollary 3 and Corollary 4, we have 



(l + r) 2 |Z0| 2 < E / \Z^\du<e 2 E . 



l/3|<3' 

To bound \d v <f)\, notice that r > R. Inequality (50) implies that 

E r I \Z p d v ip\dLudt<e 2 E 0l tp = r<j>. 
\p\< 2 Jt 2 

Then using Sobolev embedding on [ti,t 2 ] x S 2 , we obtain 

r 2 \d v ^\ 2 <\d v ^\ 2 + \^\ 2 <e 2 E . 

Therefore 

(1 + r) 2 |^| 2 < (1 + rf(\Z4>\ 2 + \d v( f,\ 2 ) < e 2 E . 
In sum, (41) follows from (43). And have proven Proposition 4. 
Remark 11. We in fact can show that 

D 1+a [z^F] T T i ^Va + n)- 1 -". 

However, it is sufficient to consider D 2a [Z^ F]?£ in order to close the boostrap argument. 
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6 Proof of the Main Theorems 

We used the foliation E T) part of which is null, in the previous argument. However, we do not have a local 
existence result with respect to the foliation E r . To solve the nonlinear equation (4) , we use the standard 
Picard iteration process. Take cj)-i(t,x) — 0. We solve the following linear wave equation recursively 

fn^ l+ i+JV($,0 n+ i) + L(50 n+ i)=F(a^ l ), 
\4> n+ i(0,x) = e<j>o(x),dt(t>n+i(0,x) = e<t>i(x). 

Now suppose the implicit constant in Proposition 4 is C\, which, according to our notation, depends 
only on R, a, t , C Q , A al3 , B af! . Set 

1 

\JC\Eq 

Then for all e < e , we have 

C\e 4 E 2 < e 2 E . 

Thus by the continuity of F{d<f> n ), we in fact have shown that the nonlinear term F satisfies 

D 2a [ZPF(d<t> n )}ll < CiE 2 e 4 (l + n)- 1 -" < E e 2 (l + n)" 1 "", V|/3| < 4, Vr 2 > n > 0. 
Therefore, Proposition 2 implies that 

E[zUn}(r) <E e 2 (l+r)- 1 - a , Vn, V|/3| < 4. 
After using Sobolev embedding on the unit sphere, Corollary 3 and Corollary 4 indicate that 

M < V^c(l + r)"'(l + |t - r + V|/3| < 2, 

I0n| < ^(l + r)- 1 . 

We also need to show that <f> n is uniformly bounded in C 2 . We first show that <fi n is bounded in C 1 . 
When r < ^i?, estimates (44) implies that 

\dZ k ^ n \ 2 < E e 2 (l + T)- 1 - a , Vfc<2. 

Here we have to point out that although there \/3i\ < l(due to the assumption \j3i\ + |/?2| < 3), the 
estimate holds for \(3i\ < 2 if we assume \(3i \ + 1/?2| < 4, see Remark 9. When < r < R, using (47) and 
Sobolev embedding on the unit sphere, we obtain the same estimates as above. For r > R, the inequality 
(52) implies that 

r-T 2 



f 2 supr-« f \dM 2 duodT < (1 + Tl )- 1 - a e 2 E , $ n = rZ^<f> ni \/3\ < 3. 

Using Sobolev embedding on S 2 x [n, r 2 ], we obtain 

|rd„Z^„| 2 < ^ 2 + r a (l + f)- 1 -^ 2 ^, V|/?| < 1. 
Recall that d u = d t - d r and |0„| 2 , |<9 t 0„| 2 < (1 + r) _1 (l + T)- 1 - a e 2 E . Wc can estimate 

\d r Z^ n \ < (l + r)-^(l + r)-^ a ^e, V|/3| < 1. 
Since = ^, we have shown that outside the cylinder {r < R} 

^(l + rr^l + Tj-H^, V|/3|<1. 



It remains to control C 2 estimates of 4> n . Outside the smaller cylinder {r < jR}, we use the equation 
(53). In fact, we can write 

2 

drr4>n+l = F{d<t> n ) + d tt (f>n+i - -d r <j> n+ i - ^(j) n +i - N($, 4> n +i) - L(d(f> n+1 ). 
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Since we have already shown that 

\dcf> n \, \d tt <j> n \ < E e 2 (l + r)"'(l + T)~^ a , 
\Hn\ < r- 2 \n 2 cf> n \ < E e 2 (l + r)"*(l + t)-*-* q , 

we conclude that 

\d rr ct> n \ < V^e(l + r)-'(l + T)-'-' a . 
Hence when r < -j, we can show that 

\8 2 (b n \ < \d rr <f> n \ + \dZ<f> n \ < v^e(l + r)-*(l + t)-i-i q 

For the case when r < -j, we rely on elliptic theory. First we have the elliptic equation for 4> n +\ 

A<j) n+1 = F{d(f> n ) + du4>n+i - N($, <f> n+1 ) - L(d<j) n+1 ). 

Since we have shown from (45) that 

\\dZ^ n \\\ <||VZ^„||V +\\Ze+^ n \\\ <E e 2 {l + r)-'-^ V|/3| < 1, 

c5 ( B iii) c5 ( B i«) c5 ( B i H ) 

we can show that the right hand side of the above elliptic equation is uniformly bounded in C?(Bi R ). 
Then Schauder estimates [6] imply that 

ll<M c ^ (Sl ^^(l + r)- 1 -. 

In particular, we have shown that 

J2 \d P M<VE~oe(l + r)-Hl + T)-^^. 

I/3|<2 

Now the classical local existence theory [8] shows that there exists a time t* > and a unique smooth 
solution <p(t,x) e C°°([0,r) x K 3 ) of equation (4). Moreover 



<j>n(t,x) -> <j>(t,x) 

in C°°([0,r) x R 3 ). Therefore 



E l^l<V^e(l + r)-'(l+T)- 1+ ' Q , V(t,aOe[0,f) 

\P\<2 



x M 3 . 



By a theorem of Hormander [7] that as long as the solution is bounded up to the second order derivatives, 
the solution exists globally. That is there exists a unique global solution <f>(t, x) <G C°°(M 3+1 ) which solves 
(4). In addition since 

4> n (t,x) -><t>{t,x), (i,i)£K 3+1 , 
<f> obeys all the estimates of <f> n obtained above. We thus finished the proof of Theorem 2. 

To prove Theorem 1, it suffices to check that the functions <E>, A/" M (<9$), N^ u {d<&) satisfy the conditions 
in Theorem 2 that L'*, W satisfy correspondingly. In fact, notice that 

\dZPN»(d<$>)\ < C{M)\d 2 Z^l V|/3|<2. 

The boundedness of d 2 Z^<^, for all |/3| < 2, follows from the equation for Z 13 ^ (similarly, express the only 
unknown term d rr ZP<& as a combination of terms with known L°° norm). Hence \dZ^N{d<^)\ + \d 2 Z^^\ 
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is bounded by a constant depending on Co and the nonlincarity N. For the other conditions, when t < to, 
we have 

|Z^"(a*)| + \Z^Af^(d^)\ < C{M,C ), V|/3| < 4. 
When t>t n , we can show 

\Z^\ < C(W,C )(1 + |ar|)- 1 -5«o(i + ( t _ | ^ | ) _,_ ) ~ 1 , 
|^A^' y ($)| < C(jV,C )(l+ M) - * 00 , |0| < 4. 

Replace a with minj^ 1 , a}. Then the functions A/" p (<9$), N^ v {d<&) satisfy the conditions in Theorem 2. 
We thus have the stability result of Theorem 1. 
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